Magnetoresistance of itinerant electrons interacting with local spins 
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Transport properties of itinerant electrons interacting with local spins are analyzed as function of 
the bandwidth, W, the exchange interaction, J, and the band filling, n, near the band edge. Nu- 
merical results have been obtained within the dynamical mean-field approximation and interpreted 
with the help of analytical treatments. If spatial correlations of the magnetic fluctuations can be ne- 
glected, and defining p(rn) as the magnetization dependent resistivity, we find that 5p/p(0) ~ n~ 4 ^ 3 
for J/W, n < 1 (weak coupling, low density limit), Sp/p(0) ~ constant, for J/W > l,n< 1 (double 
exchange, low density limit), where 5p = dp/dm 2 \ m 2_ t0 . Possible limitations from ignoring both 
localization effects and critical fluctuations are also considered. 

PACS numbers: 



I. INTRODUCTION. 

A large variety of ferromagnetic materials can be de- 
scribed in terms of a band of itinerant electrons inter- 
acting with localized spins. Among others, we can in- 
clude the manganese perovskites, La2_ :E Re :r Mn03]l|], the 
doped pyrochlores, derived from Ti2Mn^70, and 
many doped magnetic semiconductors 0, Q- Many of 
these materials show unusually large values of the mag- 
netoresistance near the Curie temperature. The simplest 
description of the electronic and magnetic properties of 
these compounds includes an itinerant band of non in- 
teracting electrons, parameterized by the bandwidth W , 
and a finite concentration of classical localized spins, 
which interact with the electrons through a local ex- 
change term, J. In the following, we assume that the 
concentration of spins is of the same order as the num- 
ber of unit cells in the lattice, and describe the model 
in terms of two dimcnsionlcss parameters, the ratio J/W 
and the number of electrons per unit cell, n. The values of 
these parameters for the pyrochlore compounds typically 
satisfy J/W,n <C 1, while the double exchange model 
appropriate for the manganites is such that J/W ^> 1, 
and n can be of order unity. The simple description used 
here is probably inadequate for the dilute magnetic semi- 
conductors, where the concentration of local moments is 
small, and, moreover, disorder effects can lead to local- 
ization of the electronic statesQ. 

It has proven useful to identify general trends in the 
transport properties of these materials. The resistance is 
typically high in the paramagnetic phase, and decreases 
as the magnetization increases. A simple parameter used 
to characterize the magnetoresistance is the dimension- 
less ratio C = [dp(m 2 ) / 'dm 2 ]/ p| m 2_>o, where p(m 2 ) is 
the magnetization dependent resistance. The resistiv- 
ity is dominated by scattering processes with wavevector 
|q| ~ 2/cf|8|) where Uf is the Fermi momentum. The 
scattering arises from magnetic fluctuations, which are 
described by the magnetic susceptibility, x(q). Assuming 



that x(q) can be approximated taking into account only 
the low momentum critical fluctuations which exist near 
the Curie temperature, one finds that C oc n~ 2 / 3 [^. ITo|. 

A different regime exists at sufficiently large values of 
the electronic density, when x(|q| = 2kp) cannot be de- 
scribed in terms of critical fluctuations. In a system with 
a large coordination number, the properties of x(q) f° r 
sufficiently large values of q can be studied using an effec- 
tive medium approach, which neglects the spatial corre- 
lations of the magnetic fluctuations. As it is well known, 
this description becomes exact when the dimensional- 
ity, or the coordination number, becomes large I fj. The 
conductivity of models of itinerant electrons interacting 
with classical spins has already been analyzed in this 
limit [13, [13J. In the following, we will study numerically 
the dependence of the parameter C defined above on the 
density of itinerant electrons, mostly in the regime when 
€f <C W, where If is the Fermi energy measured from 
the bottom of the band, assuming uncorrelated magnetic 
disorder. We will perform numerical calculations in the 
Coherent-Potential Approximation and use simplified an- 
alytical methods to understand the numerical results. We 
will find that, for J/W <C 1, the band splitting due 
to the magnetization leads to the divergent behavior, 
C oc n -4 / 3 . On the other hand, the magnetoresistance 
coefficient C remains finite close to the band edge in the 
double exchange limit J/W ^> 1, a fact easily under- 
stood in terms of a magnetization-dependent bandwidth. 
The paper is organized as follows. Section [H] presents the 
model and the perturbation limit for C. The numerical 
method is described in section ITTT1 Section llVl contains 
results for the weak coupling limit, with a simplified ana- 
lytical treatment and a discussion of possible limitations 
due to localization effects. Results for the strong coupling 
limit and a simplified explanation are found in sectionlVl 
The last section summarizes the main conclusions of our 
work and makes contact with the experimental situation, 
placing out results in the context of available compila- 
tions of magnetoresistance data0, 0, [w| . 
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II. MODEL HAMILTONIAN. 
MAGNETORESISTANCE IN THE 
PERTURB ATIVE LIMIT 

We consider free electrons coupled to core spins and 
modeled by the following Hamiltonian: 



H 



2> 



J 

-k.cr 



Ck,, 



(1) 



where a = ±1 is a spin index, is the band struc- 
ture of carriers in the unperturbed system, <Ji are Pauli 
spin operators for carriers at site i, and J is the cou- 
pling constant between carriers and core spins Si, the 
latter treated classically as unit length vectors Sf = 1. 
It should be understood from the beginning that our aim 
is the study of transport properties of the carriers, and 
not the ground state of H. Therefore, we consider the 
random distribution of core spin as given, parameterizing 
the problem. Consistently with the discussion of the in- 
troduction, we assume this probability distribution to be 
site-uncorrelated, and chosen to describe the change from 
a paramagnetic situation to a (weakly) ferromagnetic 
one. This can be achieved, for instance, with a prob- 
ability 'P(S) oc exp(h e ffS z ). Within this parametriza- 
tion, h e ff — corresponds to the paramagnetic phase, 
and h e ff ^ describes a phase where the core spins are 
polarized along the z axis. 

Following standard practice, we define the magnetore- 
sistance coefficient, C, as: 



p(0) - p(m) 
P(0) 



C to 



(2) 



where p(m) is the resistivity corresponding to an average 
polarization of core-spins m = (S z ). Notice that eq. |2] 
only makes sense to order m 2 , the only situation we will 
consider here. 

Before embarking on more complex formalisms, let us 
consider what perturbation theory (in the coupling J) 
has to say about the magnetoresistance. In the standard 
relaxation time approximation, the resistivity is given by: 



We conclude that the magnetoresistance coefficient C is 
given by 



C = l, 



(6) 



a universal number, independent of the carrier concentra- 
tion. This uninteresting results seems to preclude further 
consideration of the weak coupling limit, but, as we will 
show in the following sections, this is not the whole story. 
We will obtain that, close to band edges, this perturba- 
tive result fails, leading to an enhancement of C with 
decreasing carrier concentration. 



III. COHERENT-POTENTIAL 
APPROXIMATION 

The Hamiltonian of Eq. ^ describes non interacting 
electrons moving in a random potential. The obten- 
tion of electronic properties requires, therefore, perform- 
ing configuration averages over the random orientation 
of core spins. The Coherent-Potential Approximation 
(CPA)[l]j Hal, a procedure that we sketch here, offers 
a convenient way in order to obtain such averages (see, 
for instance, ref.^3 f° r coverage of the original litera- 
ture). We start with the first term of Eq. as our 
unperturbed Hamiltonian H D = ^ k a Ck^ a , lead- 

ing to a Green's function or resolvent whose local matrix 
elements are given by: 



,(*) = {i,a\{z-H )- l \i,a) = / de 



Vo{e) 



(7) 



where \i, a) is the electron state at site i and spin cr, and z 
is a complex energy in the upper plane. This unperturbed 
system is characterized by a reference density of states 
(per spin) T> (e) = — tt~ Im g Q (e + i0 + ), which we take 
to be the usual semielliptical function: 



2 
ttW 



(8) 
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where n is the number of carriers, e the unit charge, m,b 
the band mass, and r is the relaxation time, given to 
lowest order in the random potential by: 



T- 1 CX J 2 ((S 2 



(S) 2 ) V(e F ) 



(4) 



being T>(6f) the density of states per spin at the Fermi 
level. Upon magnetizing the core spins, only changes 
in the fluctuating potential will modify the conductiv- 
ity to order J 2 . Therefore, the relevant dependence is 
r _1 (m) ~ J 2 (l — to 2 ), leading to: 



p(0) - p(m) 
P(0) 



^CQ-t-^to) 

T-1(0) 



TO 



(5) 



notice that band edges are entirely consistent with a 3-d 
system. 

The configuration average of the Green's function for 
the entire H: g„(z) = ((i, a\(z — a)) (external 

brackets stand for configuration average) is calculated by 
the CPA from that of the unperturbed system by means 
of a self-energy E CT , in the following manner: 



9a{z) = g a {z - 



(9) 



where the self-energy, Yi a (z), itself a function of z, can 
be thought of as providing an effective medium that sub- 
stitutes the real system. This self-energy is obtained by 
imposing that the local scattering produced when we re- 
place an effective site by a real one be, on the average, 
zero 01- Adapted to our case, the self-energy obeys 
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the following equation (rotational symmetry around the 
z axis assumed): 



1 



9o(z 



= (F a (z,e)- 



J 2 sin 2 ( 

F-a(z,t 



(10) 



where cos 9 = S z , brackets stand for averages 
over the distribution of core spin orientations, 
and the auxiliary function F a (z,9) is given by 
F a (z, 9) = Ecr + aJcos9 + l/g {z — where 
a = +1(— 1) for up (down) electron spin. 

Although primarily intended for one-body properties 
(such as the density of states) , the CPA can be applied 
to transport properties. The basic idea consists in decou- 
pling the two-body correlations that appear in the Kubo 
formula into products of one-body correlators, which are 
then obtained with the CPA self-energy. Skipping further 
details, the relevant expression (see ref.^Tj) adapted to 
our case is 



P 



E^ 1 



where the conductivity for spin (<r) carriers, p a 
by: 



(11) 



is given 



Pa 1 oc 



de'f(e') Im 



1 



e F + i0+ - S CT (e F + i0+) - e' 



(12) 

where ep is the Fermi level and the function f(e') oc (1 — 
(e'/IF) 2 ) 3 / 2 accounts for matrix elements of the velocity 
operator (l7llTs| . 

The whole CPA procedure for transport properties can 
be considered as a mean-field approximation. Actually 
a double mean-field: first for the resolvent and later for 
transport (decoupling). Therefore, its main shortcom- 
ing will be fluctuation related effects (i.e. localization). 
Nevertheless, we expect the 3-d nature of our problem to 
mitigate this limitation. In fact, the CPA has been shown 
to provide a good description of 3-d disordered systems, 
being one of the very few non perturbative methods avail- 
able. In recent times, this approximation is often termed 
dynamical mean fieM|ll| , an approach invented for gen- 
uine many body problems that becomes equivalent to the 
CPA when applied to a Hamiltonian like that of eq. (TJ. 



IV. RESULTS. WEAK COUPLING LIMIT: 

J/W « 1 

We have calculated the magnetoresistance coefficient 
by direct numerical subtraction of the CPA conductiv- 
ities in the paramagnetic and ferromagnetic cases, the 
latter induced with a tiny effective field in the proba- 
bility distribution for core spins: V(S) cx exp(h e ffS z ). 
Representative results for coupling constants well below 
the bandwidth (J << W) are shown in fig. ^ where 
the magnetoresistance coefficient (C) is plotted versus 
carrier density per site (n). Notice that, in this weak 




FIG. 1: Magnetoresistance coefficient C (eq. |5J versus carrier 
density for two values of J. Inset: enlarged view of the low 
density region. 



coupling regime, the perturbative value C = 1 is closely 
approached over most of the concentration range. Yet, 
significant enhancement of C is evident for small carrier 
concentration, that is, close to the band edge. The inset 
of fig. H blows up the low density regime, where values 
of C > 100 are easily obtained. In fact, our numerical 
results are compatible with a divergent behavior of C in 
that limit. Fig. ^is representative of the small coupling 
regime J << W: no matter how weak J is chosen, we can 
always find large enhancements of C if we approach the 
band edge. In fact, the smaller the value of J, the smaller 
the range of densities required to see this enhancement, 
as shown in the inset on fig. ^ Notice that this effect 
does not formally contradict the perturbative result: for 
fixed carrier density (i.e. fixed Fermi level), C — > 1 when 
J — * 0. Nevertheless, the quantitative failure of the per- 
turbative results for a given J seems unavoidable upon 
approaching the band edges. In what follows, we will ad- 
dress the origin of this failure with a simplified analytical 
procedure. 



A. Effect of band splitting. Analytical treatment 

Our numerical results (and the analysis to follow) in- 
dicate that the nominal criterion for the validity of the 
perturbative result, J << W , valid well within the band, 
has to be replaced by J << e, where e = e F — (— W), is 
the distance of the Fermi level to the band edge. To see 
how this comes about and its effect on the magnetore- 
sistance, the original CPA equations are rather opaque 
and inconvenient. Instead, we will start with the per- 
turbative treatment of section [HI adding the effect we 
believe is at the origin of the enhancement: band split- 
ting. The first manifestation of a net polarization of core- 
spins, m, is a band splitting of the unperturbed density 
of states T> (e ± Jm). This splitting does not show up 
in the resistivity if one keeps the calculation to order J 2 , 
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FIG. 2: Schematic density of states showing the relative po- 
sitions of Fermi levels upon magnetization. 
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FIG. 3: Magnetoresistance coefficient C versus carrier density 
for two values of J illustrating the scaling behavior of eq. 1151 
in the low density limit. 



but will affect higher order terms. The situation is de- 
picted in fig. [21 where the up and down densities have 
been aligned to share a common origin, leading to dif- 
ferent apparent Fermi levels for up and down electrons: 
?! — e_i = 2m J. Transport properties depend on the den- 
sity of states at the apparent Fermi levels which, in turn, 
changes rapidly upon splitting (magnetization) close to a 
band edge. Therefore, it is not unexpected that this mere 
shift can cause a large effect on the magnetoresistance. 
We will see that this is the case. 

We assume the simplified situation described in fig. |3 
with a parabolic band T> (e) ~ vt, and apply standard 
transport theory over this already split situation. There- 
fore, we have to discriminate between up and down con- 
tributions to the conductivity: 



the magnetoresistance coefficient can be written as : 



-l n ° e 
P =>^ T <7 

— ' TUb 



(13) 



where n a , is the number of carriers for each spin species, 
of course satisfying the constraint n = ^2 a n a . The in- 
verse relaxation times are given by 

rf 1 (x J 2 [V (^)((S 2 Z ) ~ m 2 ) + V^Sl + S 2 )} 
Tl 1 ex J 2 [V (e L )((S 2 ) - m 2 ) + V (ei)(S 2 x + S 2 )] 

where the first contribution comes from transitions 
within the same spin species, and the second term ac- 
counts for spin-flip transitions (notice the different den- 
sity of states in each case). 

Going from the paramagnetic situation (Fermi level at 
e a in fig. |2J to the ferromagnetic one (Fermi levels at 
efjl); the magnetoresistance picks up contributions com- 
ing from changes in carrier spin populations n a and life- 
times T a , in addition to the standard contribution already 
described in section^] Keeping consistently terms to or- 
der m 2 in eq. I|13|) . it is only very tedious to show that 



C = 1 



19 (J 
36 IE 



(14) 



We see that it is just a matter of getting close enough to 
the band edge e < J to obtain arbitrarily large correc- 
tions to the nominal perturbative result. 

In this simplified scheme, owing to the fact that n ~ 
the magnetoresistance becomes divergent in the low 
carrier density with the following law: 



-3/2 



C - 1 



W 



(15) 



This scaling seems to be obeyed by the CPA results, as 
illustrated in fig. [3| where the result for two values of 
J « W collapse onto the same straight line. This makes 
us believe that this simple treatment: splitting plus stan- 
dard relaxation time transport, captures the essence of 
the enhancement observed in the CPA calculations. 



B. Ioffe-Regel criterion 

We have seen that, for J <C W , CPA calculations and 
the simplified treatment of Eq. (|14f> support the existence 
of a large magnetoresistance coefficient close to the band 
edge. In this section we address the validity of this result. 
The main concern comes from the fact that we have to be 
close to a band edge, where localization effects (ignored 
in our mean-field a ppr oach) are expected to play an im- 
portant role [rRIT9l l2fl||. To estimate the e nergy ra nge of 
this effect, we use the Ioffe-Regel criterion 20, 2 1 j . This 
criterion can be stated saying that quantum corrections 
to transport can be expected to be relevant when the 
mean free path diminishes to become of the order of the 
de Broglie wavelength for electrons at the Fermi level. 
This sets a characteristic distance to the band edge em 
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such that, if the Fermi level is below it, e < Iir, localiza- 
tion is expected to dominate. In our case, this criterion 
can be written as the condition: 



•ir 



(16) 



where tir is the lifetime at <Lir. Using the perturbative 
result for the lifetime, this leads to: 



(17) 



Comparing this J 4 dependence with the J 2 of Eq. Q14JI. 
we see that there is ample room for observing the en- 
hancement of C upon approaching the band edge, before 
the Ioffe-Regel limit is reached. More quantitatively, we 
can define the enhancement factor at the Ioffe-Regel limit 
by: 




FIG. 4: Magnetoresistance coefficient C (eq. |5J versus carrier 
density for J/W = oo. Inset: enlarged view of the low density 
region. 



ACi R = 



19 / J 



36 V e/fl 



with the result that: 



AC 



in. 



10" 



(18) 



(19) 



Therefore, we conclude that, for our weak coupling case 
J <C W, large increases of the magnetoresistance coeffi- 
cient C are allowed before localization effects set in. 



V. RESULTS. DOUBLE EXCHANGE LIMIT: 

J/W»l 

In this limit, the splitting between the two spin sub- 
bands is large, and we can neglect the one which is at 
high energy, ~ J. We have implemented the J = oo 
limit in the CPA equations. A similar calculation has 
been carried out before by Furukawa|l2(, but here we are 
mainly interested in results close to the band edge. In fig. 
0]we show the CPA results for the magnetoresistance co- 
efficient C as a function of carrier density. Notice a mild 
enhancement of C close to the band edge but, unlike the 
weak coupling limit, no divergent behavior is observed. 
This is more clearly seen in the inset of fig. where the 
low density region is blown up, showing a saturation of 
the magnetoresistance coefficient around C ~ 11. 

Direct use of the CPA equations to understand this 
behavior is, again, inconvenient. A simpler picture in 
order to include the effect of a finite magnetization is 
to assume that the bandwidth depends on mj^] in the 
following manner: W m — W (1 + m 2 )/2. This ansatz 
interpolates between the paramagnetic (W m — W/y/2) 
and fully ferromagnetic limit (W m — W), giving a rea- 
sonable approximation to the dependence of the Curie 
temperature on band filling |23l l24j. Assuming a semiel- 
liptical density of states and measuring energies from the 



center of the occupied band, we have: 
2 / 



V(e) 



irW„ 



€ 



(20) 



and the change in the electronic Green's function can be 
described by a self-energy, S(z, to), such that: 

g{z,m)=g [z-?,(z,m)\ (21) 

with z = e + t0 + . This self-energy is: 

1 — to 2 



E(z, to) 



2(1 + to 2 ) 



+ 



Wl 



(22) 



We now assume that the conductivity is proportional to 
n x ImS(ei?, to) -1 as done in previous sections. Then: 

<9 2 Im£(e F , m) 



C~ - 



1 



ImE(eF,0) 



di 



(23) 



In this equation we have to take into account that, if the 
number of particles is fixed, €f ~ Wy/(1 + m 2 )/2(— 1 + 
[(37r)/(4v / 2)] 2 ^ 3 ". 2 ^ 3 ) has an implicit dependence on 
Then ImS - n 2 ^W(l - 



divergent term in C as n 



lim C 

n—yQ 



< 2 )/VT 

0. Therefore, 
- constant 



i 2 , and we find no 



(24) 



This saturation of C is obtained in the full CPA calcula- 
tion (see fig. lending support to this simple picture 
of a m-dependent bandwidth in the limit J — oo. No- 
tice that, for this picture to apply, the number of carriers 
must be kept fixed upon magnetization. If the chemical 
potential were kept fixed, on the other hand, we would 
have obtained a divergent behavior: lim„^oC ~ n~ 2 / 3 '. 



VI. VALIDITY OF THE APPROXIMATIONS. 
SUMMARY. 

The analysis discussed above uses an effective medium 
approach to study the transport properties of itinerant 
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FIG. 5: Value of C for the pyrochlore compounds analyzed 
m 0> Q (dark dots) and for the magnetic materials studied 
in 9] (stars). The full line is a fit to a n -4 ^ 3 dependence. The 
broken line is the n -1 ^ 3 dependence discussed in|jj. 



electrons coupled to classical spins. As mentioned in the 
introduction, the method neglects the spatial correlations 
of the magnetic fluctuations: a potentially serious draw- 
back if the source of polarization is the spontaneous mag- 
netic order expected for the Hamiltonian of eq. The 
critical behavior of these fluctuations is particularly im- 
portant in the vicinity of the point q = 0. In addition, 
our scheme cannot take into account the localization of 
the electronic states induced by the disorder in the mag- 
netization near the band edges. Both limitations become 
important when the electronic density is lowered, as the 
relevant scattering processes are shifted towards low mo- 
menta, and the Fermi energy approaches the band edge. 
We will consider these two effects separately. 

The lack of spatial correlations is a generic feature of 
effective medium theories. These correlations modify the 
momentum dependence of the susceptibilities at low val- 
ues of the momenta. If the number of dimensions of the 
system were indeed large, these effects could be safely ig- 
nored, as the volume of the region near q = is negligible 
compared to the volume of the unit cell. We are inter- 
ested in changes in the resistivity near the Curie tem- 
perature, so that critical fluctuations at low momenta 
are always present. The length scale at which these fluc- 
tuations become important is the correlation length in 
the paramagnetic phase, £o, whereas the length scale for 
transport is A^jlCJ]. Therefore, we expect our results 
to be applicable for kp£,o ^> 1, leading to the following 



condition for the carrier density: n S> (a/Co) D > where D 
is the dimension and a the lattice size. Assuming, for 
instance, a typical value £o = 5a and D = 3, we find 
n > 1(T 2 . 

Localization becomes important when the mean free 
path associated to the disorder, I is such that kpl ~ 1. 
In the weak coupling limit, its effect has been discussed 
before, with conclusions that we briefly repeat here. Lo- 
calization is expected to be relevant below the Ioffe-Regel 
energy, em ~ J 4 /W 3 . Therefore, observing the en- 
hanced behavior C ~ rt~ 4 / 3 requires the following con- 
dition on the Fermi level: 

^3 < i F < J, (25) 

a regime easily accessible in the weak coupling limit, J <C 
W. 

In the double exchange limit, J/W ^> 1, we have 
ImE ~ y/WeF > ep for n <C 1. The mean free 
path becomes comparable to the lattice spacing at low 
densities, implying that localization effects are relevant 
in this regime. Notice, though, that the applicability 
of the previous criterion, imported from the perturba- 
tive limit, is open to question in this strong coupling 
limit (2j| . Moreover, the use of the scheme used here to 
study the phase diagram shows that the paramagnetic- 
ferromagnetic transition becomes first order for n < 
10 -1 |23|. This inclusion of this effect can change sig- 
nificantly the magnetoresistance. 

In summary, we have studied transport properties of 
electrons coupled to local, uncorrelated core spins. The 
magnetoresistance coefficient (eq. [3| has been shown to 
diverge in the weak coupling limit as C ~ n~ 4 / 3 , near the 
band edge, as a consequence of band splitting upon mag- 
netization. On the other hand, the magnetoresistance 
remains finite in the strong coupling limit (double ex- 
change model) close to the band edge. This is understood 
from a simplified picture consisting in a magnetization- 
modulated bandwidth. Possible limitations of the previ- 
ous results from the neglect of both localization effects 
and critical fluctuations have been addressed. 

Finally, we would like to consider the potential ap- 
plicability of our results to the experimental situation. 
In Fig. [5], we present available information^,^] for the 
magnetoresistance in pyrochlores and other magnetic ma- 
terials. We observe a strong dependence on the density 
for the pyrochlores family, which seems to come closer to 
the n -4 / 3 law obtained in this paper than the fit based 
on 2kp scattering alone 
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